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Abstract 

Polarized x-ray absorption near edge structure (XANES) spectra at both the 
Cu and the Ge if -edges of CuGeOa are measured and calculated relying on 
the real-space multiple-scattering formalism within a one-electron approach. 
The polarization components are resolved not only in the unit cell coordi- 
nate system (e||a, e\\b, s\\c) but also in a local frame attached to the nearest 
neighborhood of the photoabsorbing Cu atom. In that way, features which 
resist a particular theoretical description can be identified. We have found 
that it is the out-of-Cu04-plane Pz> component which defies the one-electron 
calculation based on the muffin-tin potential. For the Ge if edge XANES, the 
agreement between the theory and the experiment appears to be better for 
those polarization components which probe more compact local surroundings 
than for those which probe regions with lower atomic density. 
Paper published in Phys. Rev. B 66, 155119 (2002) and available on-line at 
http://link.aps.org/abstract/PRB/v66/el55119. 
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I. INTRODUCTION 



The main incentive for studying CuGeOa comes from the fact that it exhibits the spin- 
Peierls transition.^ Apart from that, it is a member of the very interesting family of copper 
oxides, which displays a manifold of outstanding physical properties. X-ray absorption spec- 
troscopy (XAS) is a convenient tool for investigating low-lying unoccupied electron states, 
due to its chemical and angular selectivity. This is especially important for compoimds, 
where by exploring different absorption edges, XAS is able to provide a comprehensive pic- 
ture of the local electronic structure of the material in question. Further details can be 
obtained by analyzing polarized (i.e., angular-dependent) spectra, as in that way features 
which would be obscured in unpolarized spectra can be observed. Also, comparison with 
polarized experimental x-ray absorption spectra poses a much more stringent test to the 
theory.^ 

CuGeOs has an orthorhombic crystal structure^ with a=4.80 A, 6=8.48 A, and c=2.94 A 
(space group no. 51, labeled Pbmm or D2^), therefore, its polarized spectra generated by 
dipole transitions can be decomposed into three partial spectral components.^ While all Cu 
and all Ge atoms are equivalent within the symmetry properties of the lattice, there are two 
distinct O sites in the CuGeOa crystal. Instructive pictures of the CuGeOa structure can 
be found in Fig. 1 of Ref. 5, Fig. 1 of Ref. 6 or, from a more local point of view, in Fig. 1 of 
Ref. 7. 

In the past, x-ray absorption near-edge structure (XANES) of the £\\b and e\\c polarized 
components at the Cu K edge was measured by Cruz et af" and polarized O K edge spectra 
were measured by Corradini et al.^ No measurement of the Ge K edge XANES has been 
reported so far. Cruz et aP perform a theoretical analysis of the pre-peak relying on their 
many-body cluster calculation, Corradini et aP compare the lowest in energy part of their 
spectra with densities of states around oxygen atoms as provided by many-body calculations 
of Villafiorita et al.^^ No attempts to calculate the spectra in the whole XANES range were 
made in either Ref. 8 or Ref. 9. Relying on quantitative arguments and on analogies with 
spectra of other compounds, Cruz et aP suggest a many-body interpretation of the main Cu 
K edge XANES peak splitting. 

Due to the presence of strong electron correlations in CuGeOa, a one-particle approach 
apparently cannot be fully trusted for describing the states at the bottom of the conduc- 
tion band (pre-edge region, photoclcctron energy ii^ ^ 5 cV above the threshold). It may 
be worthwhile, however, to apply the one-particle formalism based on the local density 
approximation (LDA) to the rest of the XANES region. One would be then able to comple- 
ment qualitative and empirical arguments, which have been applied for interpreting CuGeOs 
spectra in this energy domain so far, with a material-specific ab-initio approach. 

In our study, we investigate polarized XANES spectra at the Cu K edge and the Ge K 
edge, both experimentally and theoretically. By recording the spectra for several orientations 
of the crystal with respect to the x-ray polarization vector e, we are able to explore not only 
the e\\b and £\\c polarizations analyzed by Cruz et aP but also the remaining £\\a component. 
Our theoretical treatment is based on the multiple-scattering formalism in the real space. 
For a more comprehensive insight into the physical nature of the electron states probed by 
the Cu K edge XANES, we resolve the polarization components not only in the unit cell 
coordinate system but also in a local reference frame attached to the coordination polyhedron 
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of the photoabsorbing Cu atom. In that way, features which defy a particular theoretical 
description can be specifically identified. As we include also the Ge edge into our study, a 
complete view of the structure of unoccupied electron states in CuGeOs as seen from the 
copper, germanium, and oxygen^ sites emerges finally. 

II. RESULTS 
A. Experiment 

We used an experimental setup very similar to that described in a greater detail in our 
previous publications. ^'^^ The layered structure of CuGeOa allows to prepare a sample for 
x-ray transmission experiments through splitting a single-crystal slab into thin plates perpen- 
dicular to the unit cell direction a. A plate of about 7 mmx4 mmx0.03 mm was separated 
from the slab by means of usual scalpel and immediately used as the sample. The experi- 
ments were carried out at the beam lines Al and E4 (HASYLAB, DESY) equipped with an 
Si (111) two crystal monochromator. For detection of the x-ray intensities in front of and be- 
hind the sample plate, usual ionization chambers were used. The sample plate was positioned 
in a PC-controlled goniometer, allowing three perpendicular rotations. In that way, several 
polarized spectra for various orientations of the polarization vector e with respect to the 
sample were recorded. In order to separate the absorption resulting from the i^'-transitions 
exclusively, we used a subtraction of the background function f{E) = a/E^ + b/E^ + c after 
Victoreen. All the spectra were corrected for an equivalent effective sample thickness, with 
an additional matching normalization of max±5% in the remote extended x-ray absorption 
fine structure (EXAFS) region. Some of the measured curves are displayed in Fig. 1 (Cu 
edge) and Fig. 2 (Ge edge) for several representative orientations of the polarization vector 
£ with respect to the CuGe03 monocrystal. 

B. Symmetry decomposition 

Theoretical foundations needed for extracting symmetry-resolved partial spectral com- 
ponents of x-ray absorption spectra were outlined by Brouder^ in an extensive way. For an 
orthorhombic crystal, in particular, any polarized x-ray absorption spectrum can be decom- 
posed into a weighted sum of three independent partial spectral components (in the dipole 
approximation). These partial spectral components can be chosen so that they correspond 
to the polarized spectra recorded with the e vectors parallel to each of the crystal axes. 
Hence, the dipole part of any linearly polarized spectrum can be written as^'^^ 

f^D = PxSl + PySl + PzSl , (1) 

where e^, £y, and are cartesian components of the polarization vector e and p^, Py, 
and pz denote the partial spectral components (their designation refiects the corresponding 
projected density of states). Note that our Eq. (1) can be transformed into the tensor form 
presented by Brouder^ by simple algebraic manipulations. 

By inverting Eq. (1) for a triad of measured spectra, the partial spectral components in 
a unit cell reference frame can be extracted (i.e., we take x\\a, y\\b, and z\\c — see Fig. 3 
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FIG. 1. Measured XANES curves at the Cu edge of CuGeOs for several representative ori- 
entations of the polarization vector e with respect to the CuGe03 monocrystal. The projected 
components of the e vector on the crystal unit cell axes a, h and c are written in square parenthe- 
ses in the legend. 
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FIG. 2. Measured XANES curves at the Ge edge of CuGeOs for several representative polar- 
ization vectors. This drawing is analogous to Fig. 1. 
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FIG. 3. A schematic depiction of the unit cell reference frame we use for decomposing the 
polarized spectra. The unit cell is drawn in the center of the plot. Atoms belonging to eight 
CuGeOs unit cells are shown in total. The symbols 01 and 02 used in the legend refer to two 
inequivalent oxygen sites. Note that in this reference frame one has x\\a, y\\b, and z\\c, as indicated. 
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FIG. 4. (a) Experimental (lower graph) and theoretical (upper graph) Cu edge polarized 
XANES partial spectral components Px, Py, and Pz resolved in a unit cell reference frame, (b) 
Experimental (lower graph) and theoretical (upper graph) Cu K edge partial spectral components 
Px', Py', and Pz' resolved in a local reference frame attached to the nearest Cu neighborhood. The 
theoretical curves were obtained for a cluster of 167 atoms, taking into account a relaxed and 
screened core hole. 

for a schematic depiction) . The results of such a decomposition of experimental spectra are 
shown in the lower graphs of Fig. 4a for the Cu K edge and of Fig. 5 for the Gc K edge. 
We checked that the partial spectral components obtained in this way do not depend on 
the particular choice of the triad of experimental spectra which was involved in inversion of 
Eq. (1) (see also Refs. 2, 11 for a more thorough discussion). Our py and components 
of the Cu K edge XANES are in a good agreement with the corresponding £\\b and e||c 
components measured by Cruz et al.^ The remaining px component cannot be measured 
directly for a crystal cleaved parallel to the he plane and is best accessible through the kind 
of decomposition we performed here. 

By considering the quadrupole transitions as well, one might be able to separate the 
dipole and the quadrupole contributions in the pre-peaks at 8977 eV and at 11089 eV and 
even to estimate the angular character of the quadrupole parts, analogously as it was done 
for CuO."*^^ However, the pre-peak intensity is too low in this case to allow a reliable analysis 
of this kind. Therefore, we limit ourselves to exploring dipole transitions exclusively. 

In order to get a better physical insight into the nature of the some spectral features, 
one may attempt to perform the spectral decomposition in other coordinate frames as well. 
In particular, this might be especially useful in a reference frame which reflects the local 
symmetry around the photoabsorbing site. Although Eq. (1) is strictly valid in the unit 
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FIG. 5. Experimental (lower graph) and theoretical (two upper graphs) Ge K edge XANES 
partial spectral components Px, Py, and Pz resolved in a unit cell reference frame. Theoretical 

curves obtained for a cluster of 169 atoms arc shown both for a ground state potential (no core 
hole) as well as for the case when a relaxed and screened core hole is taken into account. 
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cell reference frame only, it has been demonstrated that a consistent resolving of partial 
spectral components may be possible in local reference frames as well, provided that the 
coordinate axes are chosen in a suitable way^^ One can make an a posteriori check whether 
the coordinate frame was chosen in a suitable way or not by performing the inversion of 
Eq. (1) for several independent triads of measured spectra. Only if the partial spectral 
components obtained from different sets of data coincide, the corresponding coordinate 
system is a "good" one. The physical interpretation of partial components in such a local- 
symmetry-adapted system is often more transparent than in the case of a unit cell reference 
frame and may be more suitable for comparing spectra at absorption edges of atoms of the 
same element in different compounds. 

The nearest neighborhood of Cu in CuGeOs is similar to that of other copper oxides: 
Copper is in the center of an oxygen rectangle. Four nearest oxygens at a distance of 1.93 A 
form its corners, the sides of the rectangle are 2.50 A and 2.94 A long, respectively. As can 
be seen from Fig. 3, there are actually two interpenetrating networks of CUO4 rectangles in a 
CuGeOa crystal, tilted one to another at an angle of 69°. Thus one has to take into account 
when inverting Eq. (1) that the measured spectra reflect in fact an average over Cu sites 
belonging to two different systems. A similar situation arises, e.g., in a pure copper oxide 
CuO (Ref. 11). We have found that the local partial spectral components p^', Py', and p^' 
at the Cu K edge XANES can be consistently resolved in a local reference frame where the 
z' axis is perpendicular to the CUO4 plane, the x' axis coincides with a Cu-0 bond and the 
y' axis is perpendicular to both x' and z' axes, meaning that it lies in the CUO4 plane and 
holds a 9° angle with another Cu-0 bond. A schematic drawing of the orientation of this 
local reference frame with respect to the CUO4 rectangle is shown in Fig. 6. Partial spectral 
components resolved in the local reference frame are displayed in Fig. 46. It should be noted 
that a consistent decomposition of Cu K edge XANES in CuGeOs cannot be obtained if 
one chooses the local x' and y' axes parallel to the sides of the CUO4 rectangle (the partial 
components derived from different sets of measured spectra do not coincide in that case). 
Hence, the situation is just opposite to the case of CuO, where the partial components could 
be consistently resolved in a local reference frame only if the x', y' axes were parallel to the 
sides of the CUO4 parallelogram.^^ 

We did not attempt to decompose the Ge spectra in a local reference frame as our primary 
concern is the Cu edge in relation to spectra of other copper oxides. 



C. Theory 

Theoretical Cu ii'edge and Ge if edge partial spectral components were evaluated by cal- 
culating XANES spectra with the polarization vector e being pointed to relevant directions. 
Self-consistent scattering potentials were taken over from SCF-Xa molecular calculations^^ 
for clusters of 15-24 atoms, employing an amended XASCF code of Case and Cook.^^'^^ 
The results presented in in upper graphs of Figs. 4-5 were calculated using the real-space 
multiple-scattering (RS-MS) formalism, employing clusters of 167 and 169 atoms for the 
Cu and Ge edges, respectively. These calculations were done using the RSMS code, which 
is an amended descendant of the ICXANES code^^ and is maintained by our group. A few 
more technical details about our implementation of the RS-MS technique can be found 
elsewhere. ^^'^^ In order to verify the robustness of our results, we compared the results 
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FIG. 6. The local reference frame used for resolving partial spectral components Pxi, Py', and 
Pz' in this paper. The x' and y' axes are perpendicular to each other, the z' axis is perpendicular 
to the plane of the CUO4 rectangle. 

obtained by our codes and by the selfconsistent feff8 code^^ (version 8.10). Both code 
packages rely on a bit different implementations of the one-electron multiple-scattering for- 
malism: E.g., the interstitial region is suppressed in the stage of finding the selfconsistent 
potential in the feffS code while we use a proper muffin-tin formalism and Norman muffin- 
tin sphere radii are employed in feffS while we use matching-potential radii. Despite these 
differences, both codes gave similar results for identical clusters (we used 95 atoms for Cu 
edge and 91 atoms for Ge edge spectra in this comparative study). The main difference 
consisted in up to 20% lower intensities of some low-energy {E < 25 eV) spectral peaks for 
the FEFF8 code, giving thus worse agreement with experiment than our XASCF and RSMS 
codes. We did not try to optimize feffS input parameters to identify the source of this 
difference, our goal has been just to check that the two codes do not lead to contradictory 
pictures. A comparison of XANES spectra calculated for various settings of the muffin-tin 
potential parameters can be found, e.g., in Ref. 2. 

The XANES calculations for 167/169-atoms clusters were performed both without any 
core hole and with a relaxed and screened core hole taken into account (i.e., a Is electron 
was removed from the central atom and put among the valence electrons and only thereafter 
the electronic structure of the cluster was selfconsistently calculated). Neither choice leads 
to essentially better agreement between the theory and experiment than the other one. 
Generally, the core hole increases the intensity of peaks at the low-energy part of the spectra. 
This effect is not significant at the Cu edge so we display only the results of the calculation 
which takes the core hole into consideration in Fig. 4 (the intensities of theoretical peaks at 
8993 eV for the Py component and at 8997 eV for the Pz component would be about 10% 
lower if the core hole was neglected). At the Ge edge the core hole effect is more significant 
so we display both sets of curves in Fig. 5. It appears that including the core hole actually 
"over corrects" the ground state results in this case and consequently does not give rise to 
distinctly better agreement between theory and experiment for the Ge edge XANES. 

We made analogous calculations for non-selfconsistent potentials constructed according 
to the Mattheiss prescription (superposition of potentials and charge densities of isolated 
atoms), too. As could be anticipated, selfconsistent potentials give rise to theoretical spectra 



9 



which agree better with experiment than in the case of non-selfconsistent potentials, however, 
the improvement is not an essential one. All distinct spectral features that can be observed 
at the theoretical curves in Figs. 4-5 are present also in corresponding spectra obtained 
for non-selfconsistent potentials, with peak positions and intensities differing no more than 
Ri20% (we do not show those curves for brevity). Thus use of self-consistent potentials is 
not crucial for our study. 

The alignment in energy of the theoretical curves shown in Figs. 4-5 was set so that 
the best overall agreement between theory and experiment is achieved. All theoretical 
curves were convoluted with a Lorentzian function of the width w given by the ansatz 
w — Wc + 0.03 X E, where the constant part Wc accounts for the core hole hfetime^° and the 
energy-dependent part mimics the finite lifetime of a photoelectron with an energy E. The 
factor 0.03 was chosen just by convenience — we did not attempt to achieve the best possible 
agreement between theory and experiment by optimizing the smearing function. 

III. DISCUSSION 

Let us concentrate first on the Cu K edge in the unit cell reference frame, where a 
partial comparison with the work of Cruz et af is possible. As it can be seen from Fig. 
4a, our calculation describes the gross shape of the three components p^, and their 

polarization splitting. However, only the pz component, probing states parallel to the c 
axis, is reproduced by the theory as concerns the number of features (shoulders at 8986 eV 
and 8993 eV and peaks at 8997 eV, 9000 eV, and 9008 eV at the experimental curve), 
their positions (within 1-2 eV's), and relative intensities. The reason why the two peaks of 
the main doublet in p^ are not clearly separated in the experiment (either ours or that of 
Cruz et afi) is not clear. It may be a consequence of experimental energy resolution (most 
probably arising from the so called thickness effect — absorption maxima are suppressed and 
consequently smeared for thick samples) or it may stem from a smearing caused by many- 
body effects, not accounted for by our theory. The other two experimental components p^, 
Py seem to include an ingredient which defies our theoretical description substantially more 
than is the case of the p^ component. 

Cruz et af interpret the Cu K edge XANES of CuGeOa in this energy range by semi- 
quantitative arguments, relying on comparison with Cu spectra of La2Cu04 and CuO. They 
suggest an essentially many-body interpretation of dominant spectral features: The double 
structure of the main peaks at 8987 eV and 8993 eV for the Py component and at 8997 eV 
and 9000 eV for the p^ component (see Fig. 4a) is ascribed to transitions to well-screened 
final states followed by satellite transitions to poorly screened states. The success of our 
one-electron LDA-type computation to describe at least the Pz component correctly seems 
to contradict this point of view. Of course, no definite conclusion regarding the one-electron 
or many-body nature of the main Pz peak splitting can be made until a proper many-body 
calculation is done. 

As noted above, the experimental Px and Py components in Fig. 4a display more serious 
deviations from the one-electron theory than the Pz component. In order to identify the 
source of these deviations, it might be useful to explore the partial spectral components in 
the local reference frame. Indeed, if we compare the theoretical and experimental curves in 
Fig. 4b, it emerges clearly that the primary source of discrepancy between the theory and the 
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experiment lies mainly outside the CUO4 plane. The peak positions and relative intensities 
of the Px' and py/ components are reproduced by theory with a similar accuracy as was the 
case with the p^ component in the cell reference frame. The absence in the experiment of a 
well-defined doublet structure of the main peak at 8997-9004 eV is probably caused by the 
same mechanism as a similar situation with the cell Pz component mentioned above. The 
local Pz' component is described by the theory significantly worse than the local px' and pyi 
components, especially as concerns the fine structure of the peak at 8988-8992 eV in Fig. 
46. 

Having this in mind, one can re-assess the situation in the unit cell reference frame: 
As the c axis runs parallel to the CUO4 plane, the Pz component probes solely states lying 
within this plane and, hence, can be correctly characterized by the theory (cf. Fig. 3). 
On the contrary, the Px and py components probe states lying within the CUO4 plane as 
well as states reaching perpendicular to it. It is the perpendicular-to-Cu04 admixture into 
these components which spoils the agreement between theory and experiment. Without this 
mixing, one could have expected for Fig. 4a a similar agreement between the theory and 
experiment as can be seen in Fig. 46 for the Px' and Py> components. 

The failure of the theory to describe the pz' component in the local reference frame is 
not a total one: The theory accounts for the general trends of the Pz' spectral curve such as 
existence of a well-separated maximum at 8988-8992 eV followed by a broad structure around 
9008 eV. The most notable deficiency of the theory for this polarization is not reproducing the 
double-peak structure at 8988 eV and 8992 eV. This failure then transforms into deficiencies 
in describing the unit cell px and Py components in this energy region. Secondly, the fine 
structure of the second broad feature (positions and relative intensities of sub-maxima at 
9003 eV and 9008 eV) is not correctly reproduced by our calculation. 

At this point, it is difficult to identify the reason of this deficiency of the theory more 
specifically. One could make use of the analogy with La2Cu04, where the theory also 
describes the states within the Cu-0 layers fairly accurately but fails in describing XANES 
components perpendicular to that layer. ^-"^'^^ Guo et af^ attribute this failure to the presence 
of shake-up satellite peaks in La2Cu04 spectra and present plausible arguments for their 
view. However, the situation may be different for CuGeOa: Among others, the failure of 
the theory for the out-of-plane XANES components appears to be much more dramatic for 
La2Cu04 than for CuGeOs (a whole peak missing in La2Cu04 spectrum vers, an inverted 
double-peak intensity in CuGcOs spectrum). Since the calculations of La2Cu04 spcctra^^'^^ 
as well as our RS-MS calculation of CuGeOs XANES rely on the muffin-tin approximation, 
one cannot rule out full-potential effects as the main culprits. Intuitively one would expect 
non-muffin-tin corrections to be more important for the more loose out-of-plane states than 
for the states in the more closely-packed plane of CUO4 rectangles, so this explanation would 
be in accord with our results. Anyway this question seems to remain open. 

As can be seen from Fig. 4, the pre-peak at the Gu K edge appears to be reproduced 
at the correct position even if quadrupole transitions are neglected in the calculation. We 
verified that if the quadrupole transitions are included, they are completely negligible in the 
whole energy range except for the pre-edge region around 8977 eV, where they contribute 
with approximately with the same intensity as the dipole transitions. This might serve as an 
indication of a mixed dipole-quadrupole nature of the pre-peak at the Gu ii'edge of GuGeOa, 
however, one should be quite cautious about this issue as the successful reproduction of the 
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pre-peak position by our calculation may be just a fortuitous coincidence: The pre-edge 
XANES intensity crucially depends on the energy position of the onset of unoccupied states 
and on the structure of electronic states very close to this onset. This kind of information 
cannot be reliably obtained via an LDA-type calculation for a strongly-correlated material 
such as CuGeOa (see Refs. 23, 24 for LDA calculations of electronic structure of CuGeOa 
and Refs. 6, 7 for calculations going beyond LDA). Note that, e.g., in the case of CuO the 
agreement between theory and experiment in the pre-edge region is seemingly better for 
a non-self consistent potential than for a self-consistent one.^*^ So we think that we cannot 
draw really reliable conclusions regarding the pre-edge region and, consequently, do not even 
display the corresponding "quadrupole-included" curves here. 

The Ge K edge is decomposed in the unit cell reference frame in Fig. 5. As noted 
already in Sec. II C, the best agreement between theory and experiment would be attained 
for curves lying somewhat between the data obtained when the core hole is present and 
when it is neglected. Both calculations account for the strong polarization (i.e., angular) 
dependence of the experimental curves fairly well. However, not all components are described 
by the theory with the same accuracy. The main peak of the component is too broad in 
the calculated spectra and neither the fine oscillations seen in the experimental Px curve are 
reproduced by the theory appropriately (though including the core hole improves this a bit). 
The main doublet in the Py component is reproduced, however, the separation of its peaks 
is 5 eV in the experiment and 6 eV in the theory. At the p^ component, the calculation 
without the core hole reproduces the 11101 eV and 11110 eV experimental peaks both in 
positions and relative intensities, just the 11105 eV shoulder is missing in the calculated 
curve. If the core hole is included, the low-energy shoulder at 11099 eV of the theoretical 
Pz curve is drastically increased, worsening thus the agreement with experiment. It seems 
therefore that a proper description of the core hole effect at the Ge edge of CuGeOs has to 
go beyond the static relaxed and screened model. We did not try to employ other simple — 
less frequently used — prescriptions for core hole treatment (unrelaxed and/or unscreened 
hole) as they have not worked well with other non-metallic systems. 

Overall it appears that the best agreement between theory and experiment has been 
achieved for the e||c component p^, next comes the e||6 curve and, finally, worst agreement 
occurs for the £\\a component p^. Although we are unable to find unambiguously what is 
the reason for this "hierarchical" behavior, we believe that the degree of packing of the 
neighborhood of Ge atoms in various directions may be a factor. Namely, the Ge atoms 
are nearly tetrahedrally coordinated by four O atoms at 1.74-1.79 A and the a, b, and c 
directions hold different angles with those oxygens. The a direction avoids the tetrahedron 
corners as much as possible, holding a 54° angle with all the Ge-0 bonds. The minimum 
angle between the b and c directions and any of the Ge-0 bonds is 36°, on the other hand. 
Moreover, the c direction points exactly to the Ge atoms which form the second coordination 
shell around germanium at 2.94 A. Hence, a rule of thumb emerges that the more open is 
the local geometry probed by a particular spectral component p^, Py, or p^, the worse 
agreement between theory and experiment at the Ge K edge arises. A possible explanation 
of this correlation might rest with the muffin-tin approximation we employ, as non muffin-tin 
effects are expected to be more significant at open geometries than at close ones. 

Experimental XANES spectra at the O K edge of CuGeOs were presented by Corradini 
et aP and by Agui et aP^ for energies up to 25 eV above the absorption threshold. The 
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interpretation of O edge spectra is complicated by the fact that there are two inequivalent 
oxygen sites in CuGeOs, meaning that one sees a superposition of two edges in the exper- 
iment actuaUy We attempted to reproduce the polarization dependence of the O K edge 
spectra measured by Corradini et aP by a RS-MS calculation for the same self-consistent 
potential which we employed for calculating the curves displayed in Figs. 4-5. The energy 
shift between the edges of the two inequivalent oxygens was set by hand, as there are several 
factors affecting its value which cannot be described with a sufficient accuracy by an LDA- 
type calculation (a brief discussion of those effects can be found, e.g., in Ref. 9). We have 
found that there is no suitable value of such a shift which would lead to the reproduction 
of the two distinct polarization-dependent peaks within the first 3 eV above the O K edge 
threshold, as observed in the experiment.^ On the other hand, the relatively broad main 
peak at 5-10 eV above the threshold does not exhibit any distinct angular-dependent fea- 
tures at all and, consequently, can be at least roughly reproduced for several different values 
of the energy shift between the two O edges, offering thus no clue for its proper value. We 
do not display the results here for brevity, concluding that our calculation cannot be used as 
a guide for interpreting the polarization dependence of the O K edge absorption spectrum of 
CuGeOa- Strong electron correlations, unaccounted for by the LDA framework, seem to be 
the most probable explanation for the failure of our calculations to reproduce the polarized 
O K edge XANES close to the threshold. Note that for that reason (unrehabihty of LDA at 
the x-ray absorption threshold in CuGeOa) we do not discuss the pre-peaks at the Cu and 
Ge edges in detail either. On the other hand, those parts of the spectra which correspond 
to photoelectron energies higher than E ^ 5 eV can be at least partially described by a 
one-electron LDA- type approach, as demonstrated above (Figs. 4-5). 

The fact that a relatively simple theory appears to be able to reproduce gross trends in 
XANES of an insulating, probably charge-transfer oxide is quite surprising on its own. It 
seems that strong electron correlations in this type of compound need not be so crucial at 
photoelectron energies above the pre- edge, contrary to what has been generally anticipated 
so far. 

IV. CONCLUSIONS 

The main features of polarized Cu K edge and Ge K edge XANES of CuGeOs can be 
described fairly well within the one-electron framework. Contrary to earher suggestions,^ the 
main peak splitting for the s||c polarization at the Cu edge can be correctly reproduced 
within the one-electron LDA-type approach (no need for superpositioning of well-screened 
and poorly-screened states). The source of the remaining discrepancies between theory and 
experiment at the Cu Kedge XANES can be better identified and understood if the polarized 
spectra are decomposed in the local reference frame attached to the nearest Cu neighborhood 
rather than in the unit cell reference frame. We have found that it is the out-of-Cu04-plane 
Pz' component which defies the one-electron RS-MS calculation for a muffin-tin potential. 
For the Ge K edge XANES, the agreement between the theory and the experiment appears 
to be better for those polarization components which probe more compact local surroundings 
than for those which are directed to the interstitial parts of the crystal. Including the core 
hole affects teh Ge edge theoretical spectra significantly more than at the Cu edge and, to 
a certain degree, leads to a worse agreement with Ge edge experiment in comparison with 
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the case when the core hole is neglected. 



ACKNO WLED CEMENTS 

The theoretical part of this work was supported by grant 202/02/0841 of the Grant 
Agency of the Czech Republic. The experimental part was supported by Hamburger Syn- 
chrotronstrahlungslabor HASYLAB Project No. 11-98-058. The use of the CRYSTIN struc- 
tural database^^ was financed by grant 203/02/0436 of the Grant Agency of the Czech 
Republic. 



14 



REFERENCES 



^M. Hase, I. Terasaki, and K. Uchinokura, Phys. Rev. Lett. 70, 3651 (1993). 
^O. Sipr, A. Simunek, S. Bocharov, Th. Kirchner, and G. Drager, Phys. Rev. B 60, 14115 
(1999). 

^ S. Brauninger, U. Schwarz, M. Hanfland, T. Zhou, R. K. Kremer, and K. Syassen, Phys. 

Rev. B 56, R11357 (1997). 
^C. Brouder, J. Phys.: Condens. Matter 2, 701 (1990). 

^ F. Parmigiani, L. Sangaletti, A. Goldoni, U. del Pennine, C. Kim, Z.-X. Shen, A. 

Revcolcvschi, and G. Dhalenne, Phys. Rev. B 55, 1459 (1997). 
^ S. Atzkern, M. Knupfer, M. S. Golden, J. Fink, A. Hiibsch, C. Waidachcr, K. W. Becker, 

W. von der Linden, M. Weiden, and C. Geibel, Phys. Rev. B 64, 075112 (2001). 
^S. Zagoulaev and L L Tupitsyn, Phys. Rev. B 55, 13528 (1997). 

^D. Z. Cruz, M. Abbate, H. Tolentino, P. J. Schiling, E. Morikawa, A. Fujimori, and J. 

Akimitsu, Phys. Rev. B 59, 12450 (1999). 
^V. Corradini, A. Goldoni, F. Parmigiani, G. Kim, A. Revcolevschi, L. Sangaletti, and U. 

del Pennino, Surf. Sci. 420, 142 (1999). 
^°A. Villafiorita, F. Manghi, F. Parmigiani, and G. Galandra, Solid State Gommun. 104, 

301 (1997). 

11 S. Bocharov, Th. Kirchner, G. Drager, O. Sipr, and A. Simunek, Phys. Rev. B 63, 045104 

(2001). 

12 k. H. Johnsson, Adv. Quantum Ghem. 7 143 (1973). 

1^ M. Gook and D. A. Gase, computer code xascf. Quantum Ghemistry Program Exchange, 

Indiana University, Bloomington, Indiana (1980). 
1^0. Sipr, F. Rocca, and G. Dalba, J. Synchrotron Rad. 6, 770 (1999). 
i^D. D. Vvcdcnsky, in Unoccupied Electronic States, edited by J. G. Fuggle and J. E. In- 

glesfield (Berlin, Springer, 1992), p. 139. 
16 D. D. Vvedensky, D. K. Saldin, J. B. Pendry, Gomput. Phys. Gommun. 40, 421 (1986). 
1^0. Sipr, computer code RSMS, Institute of Physics AS GR, Prague (1996-1999). 
1^0. Sipr and A. Simunek, J. Phys.: Gondens. Matter 13, 8519 (2001). 
i^A. L. Ankudinov, B. Ravel, J. J. Rehr, and S. D. Gonradson, Phys. Rev. B 58, 7565 

(1998). 

2° F. Al Shamma, M. Abbate, and J. G. Fuggle, in Unoccupied Electron States, edited by J. 
G. Fuggle and J. E. Inglesfield (Springer, Berlin, 1992), p. 347. 

21 J. Guo, D. E. Ellis, G. L. Goodman, E. E. Alp, L. Soderholm, and G. K. Shenoy, Phys. 
Rev. B 41, 82 (1990). 

22 G. Li, M. Pompa, A. Gongiu-Gastellano, S. Delia Longa, and A. Bianconi, Physica G 175, 
369 (1991). 

23 L. F. Mattheiss, Phys. Rev. B 49, 14050 (1994). 

2^ Z. S. Popovic, F. R. Vukajlovic, and Z. V. Sljivancanin, J. Phys.: Gondens. Matter 7, 

4549 (1995). 

2^0. Sipr, P. Machek, A. Simunek, J. Vackar, and J. Horak, Phys. Rev. B56, 13151 (1997). 
26 a. Agui, J.-H. Guo, G. Sathe, J. Nordgren, M. Hidaka, and 1. Yamada, Solid State Gom- 
mun. 118, 619 (2001). 

2'''G. Bergerhoff, R. Hundt, R. Sievers, and I. D. Brown, J. Ghem. Inform. Gomput. Sci. 23, 
66 (1983). 



15 



